Spatial correlations in the dynamics of glassforming liquids: Experimental 
determination of their temperature dependence 
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We use recently introduced three-point dynamic susceptibilities to obtain an experimental de- 
termination of the temperature evolution of the number of molecules, A'corr, that are dynamically 
correlated during the structural relaxation of supercooled liquids. We first discuss in detail the 
physical content of three-point functions that relate the sensitivity of the averaged two-time dynam- 
ics to external control parameters (such as temperature or density), as well as their connection to 
the more standard four-point dynamic susceptibility associated with dynamical heterogeneities. We 
then demonstrate that these functions can be experimentally determined with a good precision. We 
gather available data to obtain the temperature dependence of A'^corr for a large number of super- 
cooled liquids over a wide range of relaxation timescales from the glass transition up to the onset 
of slow dynamics. We find that A'corr systematically grows when approaching the glass transition. 
It does so in a modest manner close to the glass transition, which is consistent with an activation- 
based picture of the dynamics in glassforming materials. For higher temperatures, there appears 
to be a regime where A'corr behaves as a power-law of the relaxation time. Finally, we find that 
the dynamic response to density, while being smaller than the dynamic response to temperature, 
behaves similarly, in agreement with theoretical expectations. 

PACS numbers: 64.70.Pf, 05.20.Jj 



I. INTRODUCTION 

One of the puzzling questions concerning the glass 
transition of liquids and polymers is the fact that the dra- 
matic temperature dependence characterizing the slow- 
ing down of relaxation and viscous flow occurs with no 
significant variation of the typical structural length asso- 
ciated with static pair correlations (see however For 
instance, the static structure factor of a liquid barely 
changes from the high-temperature "ordinary" liquid 
down to the glass, while the main (a) relaxation time 
may rise by up to fifteen orders of magnitude. Yet, on 
intuitive ground, it is tempting to ascribe the strong tem- 
perature dependence of the dynamics and the ubiquity 
of the phenomenon, irrespective of molecular details, to 
a collective or cooperative behavior characterized by a 
lengthscale that grows as one approaches the glass tran- 
sition. This is indeed the core of most theories of the 
glass transition 0, , starting with the picture of "coop- 
eratively rearranging regions" put forward by Adam and 
Gibbs in the mid-sixties [3]. 
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From quite general arguments, one expects that a 
seemingly diverging, or at least long, timescale must be 
accompanied by a seemingly diverging, or at least large, 
lengthscale [5|. Tracking such a lengthscale has been a 
long-standing goal of studies on glassforming materials. 
An exciting development in this direction has been the 
observation and description of the heterogeneous nature 
of the dynamics in viscous liquids and polymers 0, 0, ■ 
Experiments carried out at temperatures near the glass 
transition and computer simulations on model liquids 
have shown that the dynamics of the molecules is indeed 
spatially correlated, with regions in space being charac- 
terized by mobilities or relaxation times substantially dif- 
ferent from the average for a timescale at least equal to 
the main relaxation time. These dynamic heterogeneities 
come with (at least) one lengthscale that is typical of the 
correlations in space of the dynamics. 

Measuring the lengthscale(s) associated with the spa- 
tial correlations of the dynamics is however a difficult 
task. Experimentally, the most direct measurements 
have been performed by means of solid-state multidimen- 
sional nuclear magnetic resonance, providing estimates 
of a few nanometers (between 5 and 10 molecular diame- 
ters) [1, 0, H, [13, [m ■ Unfortunately, these measurements 
can only be done in a narrow temperature range close 
to the glass transition, thereby precluding any study of 
the temperature dependence. Other types of dynamic 
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lengthscales have been extracted from experimental data, 
mostly from crossover phenomena (e.g., between Fickian 
diffusion and non-Fickian translational motion ) or by 
some sort of dimensional analysis (e.g., from the product 
of a diffusion coefhcient by a relaxation time or by the 
viscosity [l^l; see also [31 and the discussion below). 

Computer simulations give the opportunity to study 
the properties of the dynamic heterogeneities in more 
detail, with the restriction however that only the mod- 
erately viscous regime, many orders of magnitude away 
from the glass transition, is accessible. Sizes character- 
istic of one form or another of space-time correlations 
(strings [isl . microstrings [l^, compact "democratic" 
clusters etc.) have been determined. On top of those 
studies, a more generic means of extracting a length- 
scale associated with the dynamic heterogeneities, in- 
spired by the mean-field theory of spin glasses has 
been proposed and quite thoroughly investigated: the 
spontaneous fluctuations around the average dynamics 
give rise to a four-point time-dependent correlation func- 
tion, X4(t), called dynamic susceptibility [l9j : its magni- 
tude defines a "correlation volume" that follows a non- 
monotonic variation with time t and displays a maxi- 
mum for a timescalc of the order of the a-relaxation 
time [23,[Illl2, 23, 24|]. This correlation volume shows, 
in numerical simulations, a significant increase as temper- 
ature is decreased. At the present time, however, the dy- 
namic susceptibility X4(i) cannot be experimentally mea- 
sured in glassforming liquids and polymers, which pre- 
vents any overlap between simulations and experiments. 
This is because measuring Xii^) iii principle amounts 
to resolving the dynamic behavior of molecules in both 
space and time, which has proven possible up to now 
only in systems such as colloidal suspensions [H, HI], 
granular systems (27l . [28j . and foams [1^. However, a 
promising path is to measure the non-linear response of 
glass forming liquids, for example their non-linear dielec- 
tric constant, which contains an information similar to 

xiit) MM- 

In spite of the progress made in the last decade, as 
briefly surveyed above, many important questions remain 
unanswered. To list a few: is there a unique lengthscale 
associated with the dynamic heterogeneities? Is the same 
lengthscale involved in the slowing down of the /? and a- 
relaxation? What is the temperature dependence of the 
dynamic lcngthscale(s) over the whole range from the or- 
dinary liquid phase to the glass? Is there a difference of 
behavior between the "strong" and the "fragile" glass- 
formers? What kind of underlying static length, if any, 
may drive the growing dynamic length(s) in liquids and 
polymers ? 

The purpose of the present paper is not to answer, nor 
even to address, all of the above questions. We focus on 
the experimental determination of the temperature de- 
pendence of a length (more precisely, a volume or a typ- 
ical number of molecules, iVcorr) that characterizes the 
space-time correlations in glassforming liquids. We do so 
by building upon the recent proposal made by some of 



us [3^ . m, [131 . We study the experimentally accessible 
three-point time-dependent correlation functions defined 
as the response of a two-point time-dependent correla- 
tor to a change of an external control parameter, such as 
temperature, pressure, or density. In Ref. [3^ it has been 
shown that the three-point correlation functions (that we 
will equivalently call "three-point dynamic susceptibili- 
ties") can be used to provide a lower bound to the dy- 
namic susceptibility X4(i)- More recent work (33l.[33.[3q 
has in fact established that these three-point susceptibil- 
ities contain bona fide information on the heterogeneous 
nature of the dynamics. In the present article, (i) we an- 
alyze dynamic data on a variety of glassforming liquids to 
obtain three-point dynamic susceptibilities; (ii) we stress 
the arguments (and the conditions) that allow one to 
extract a number of dynamically correlated molecules 
-^cori (i) (or a correlation volume); (iii) we study the tem- 
perature dependence of A^corr for a timescale of the or- 
der of the a-relaxation time Tq, and discuss its possible 
connection with the "fragility" (i.e., the degree of depar- 
ture of Tq, from an Arrhenius temperature dependence) 
of the glassformers; (iv) we compare the relative contri- 
butions of temperature and density fluctuations to dy- 
namical correlations, (v) finally, we discuss the connec- 
tion between the extent of spatial dynamic correlations 
and some aspects of phenomenological descriptions of the 
glass transition, such as the degree of "cooperativity" of 
the relaxation processes and the nonexponential charac- 
ter of the relaxation functions. 



II. EXPERIMENTAL ACCESS TO DYNAMIC 
SPATIAL CORRELATIONS 

We have already mentioned that the natural quantity 
characterizing the spatial fluctuations of the dynamics, 
namely the dynamic susceptibility X4(0i has not so far 
been experimentally accessible in glassforming liquids. 
However, three-point correlation functions can be ob- 
tained as the response of a two-point correlator, com- 
monly determined in experimental studies of the dynam- 
ics in glassformers, to a change in an external control pa- 
rameter. Many aspects of the problem have been studied 
in great detail in Refs. [32., jJ3„^]. Here, we mainly focus 
on the interpretation of the three-point dynamic suscep- 
tibilities and on the extraction from it of a lengthscale, 
or to the least of a volume or a number of molecules. 



A. Preliminary: A phenomenological description 

Consider a dynamic process that proceeds through 
thermal activation and is characterized by a typical time, 
T, which follows an Arrhenius temperature dependence, 
T oc exp(A/fcsT), where A is the activation energy bar- 
rier, ks the Boltzmann's constant, and T the tempera- 
ture. An example is provided by the diffusion of vacan- 
cies in a solid. The activation barrier may be due to very 
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local interactions, but in any case it is crossed because 
of energy fluctuations coming from the environment. A 
simple reasoning that allows one to derive a "correla- 
tion volume" associated with the activated dynamics is 
as follows. One just asks what is the volume V* that is 
necessary to generate typical energy fluctuations of the 
order of the energy barrier A? Typical energy fluctua- 
tions in a volume V are given by the square root of the 
variance of the total energy E in the NVT ensemble, 
i^SE'^^y = pksT'^CyV, where c„ is the specific heat per 
particle in kg units and p is the mean density. Equating 
A with {SE'^)y gives immediately 



pV* 



A2 



(1) 



If Cy does not vary much or even goes to zero as T goes 
to zero (as in a solid), then V* increases at least like 
T^^ as temperature is lowered and ultimately diverges 
at T = 0. This result looks intriguing, but the "correla- 
tion volume" so defined appears more as the result of an 
exercise in dimensional analysis than motivated by de- 
scribing spatial correlations in the dynamics. Actually, 
the nature of what is exactly correlated in the volume V* 
is unclear (if not misleading). 

Despite the obvious weakness of the approach, one may 
try to extend it to describe the slow, presumably ther- 
mally activated, dynamics in viscous liquids. There, how- 
ever, the activation energy is not constant. The temper- 
ature dependence of the main (a) relaxation time r„ is 
"super- Arrhenius" and one must choose an operational 
procedure to define an effective, temperature-dependent 
activation energy. One way is to define it through the re- 
lation Ta oc exp(A(r)/A:BT), another possibility being to 
take the derivative of logTa with respect to 1/T, namely 
A(r) — fe-B ^g°f/T)'' • These two quantities are very dif- 
ferent in fragile glassformers with a pronounced super- 
Arrhcnius behavior, since A(T) = A{T)-dA{T)/d\ogT. 
Anyhow, arbitrarily choosing at this stage the second def- 
inition and repeating the above argument leads to 



pV* = 

Cy \d\OgT ^ y 



(2) 



which reduces to Eq. ([T]) when the dependence of Tq is 
simply Arrhenius. This again results in a volume V* ov 
a number of molecules N* = pV* which increases faster 
than as T decreases and approaches the glass tran- 
sition. [This "correlation volume" is not to be confused 
with the "activation volume" sometimes introduced to 
characterize the pressure dependence of the a-relaxation 
and defined as Td\ogTa/dP\T-] A variant of the above 
reasoning was in fact proposed and further investigated 
by Donth [3, with the replacement of (l/cy) by 
|A(l/ci,)| Acy/cy, where Ac^ is the jump in the heat 
capacity at the glass transition. In Donth's interpreta- 
tion, V* (or N*) is taken, with no justification, as the size 
of the "cooperatively rearranging regions" introduced by 
Adam and Gibbs 



The same line of reasoning can be followed in the NPT 
ensemble with the activation energy replaced by an acti- 
vation enthalpy, leading to a correlation volume 



^ c„ \dlogT 



(3) 



It seems clear that for this approach to have any true 
physical meaning, it needs to be put on a much firmer 
basis and insight should be provided about the nature of 
the dynamic correlation. This is what we address below. 



B. Multi-point dynamic correlation functions and 
susceptibilities 

As was shown in Ref. [s^j and further detailed in 
Refs. [H, H^l, multi-point space-time correlations in 
glassformers may be experimentally accessible by study- 
ing the response of the dynamics to a change in a control 
parameter, i.e., by inducing fluctuations in some conju- 
gate quantity such as energy or density and monitoring 
its effect on the dynamics of the system. To provide 
a self-contained presentation, we briefly recall here the 
main arguments. 

Standard experimental probes of the dynamics in liq- 
uids give access to the time-dependent auto-correlation 
function of the spontaneous fluctuations of some ob- 
servable 0{t), F{t) = {SOiO)SO{t)), where SO{t) = 
0{t) — (O) represents the instantaneous value of the fluc- 
tuations of 0{t) from its ensemble average (O). One 
can think of F{t) as being the average of a two-point 
quantity, C(0,<) = SO{0)6O{t), characterizing the dy- 
namics. This quantity fluctuates around its mean value, 
and information on the amplitude of those fluctuations 
is provided by the variance (SC{0, i)^), where SC{Q, t) = 
C(0, t) — F{t). This allows to interpret X4(0 ^ precise 
measure of dynamic fluctuations or heterogeneities, since 
N (^SC{0,t)^) = Xiit), where TV is the total number of 
particles in the system. The associated spatial correla- 
tions show up more clearly when considering a "local" 
probe of the dynamics, like for instance an orientational 
correlation function measured by dielectric or light scat- 
tering experiments, which can be expressed as 



C(0,i) = -i J dVc(r;0,t), 



(4) 



rely 



where V is the volume of the sample and c(r; 0, t) char- 
acterizes the dynamics between times and t around 
point r. For example, in the above mentioned case 
of orientational correlations, c(r; 0, t) cx ^ X^fj^i ^(^ ~ 
ri)Y{^li{0))Y{nj{t)), where fli denotes the angles de- 
scribing the orientation of molecule i, ri(0) is the position 
of that molecule at time 0, and Y{n) is some appropri- 
ate rotation matrix element; the "locality" of the probe 
comes from the fact that it is dominated by the self-term 
involving the same molecule at times and t or by the 
contribution coming from neighboring molecules. 
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The dynamic susceptibility can thus be rewritten as 
Xiit) =p [ d\ (<5c(0; 0, t)<5c(r; 0, t)) , (5) 



where the statistical translational invariance of the liq- 
uid has been taken into account and p = N/V denotes 
the mean density. The above equation shows that Xii^) 
measures the extent of spatial correlation between dy- 
namical events between times and t at different points 
of the system, i.e., the spatial extent of the dynamic het- 
erogeneities over a time span t. Thanks to a number of re- 
cent theoretical and numerical studies [20., 21, 22. 23, 2^. 
this point is now well documented. 

Unfortunately. X4(^) has not so far been measurable 
in glassforming liquids and polymers. Furthermore, its 
physical interpretation is obscured by the fact that X4 
turns out to depend strongly on microscopic dynamics 
and thermodynamical ensemble (see [13, [3 ) . A related 
and experimentally feasible route to space-time correla- 
tions in glassformers consists in monitoring the response 
of the average dynamics to an infinitesimal change, say, 
of temperature. In the case of Newtonian dynamics, such 
as the molecules in a liquid, one can show that this re- 
sponse is given through a fluctuation-dissipation relation 
by a cross-correlation function involving the spontaneous 
fluctuations of the dynamics and those of the energy [s^ . 
For simplicity, let us first focus on the NVT ensemble. 
In this case one finds 1321: 



dF{t) 



dT 



N 



N,V 



{5E{Q)5C{Q,t))^y^, 



(6) 

where 5E{{)) = E{0) — (E) is the fluctuation of the en- 
ergy per molecule at time t = 0. The subscript NVT 
means that the statistical average is performed in the 
NVT ensemble. Expressing the fluctuations in terms of 
local quantities (which is always possible [H, [13, [13| ) as 
before leads to 



d^r{6e{O;0)6c{r;0,t)) 



NVT ' 



(7) 



where we have introduced e(r; t) as the energy density 
per molecule (hence, NE{t) = pj (Pr e{r;t)). Note that 
the correlation functions in Eqs. ^\7} are expected to be 
negative since an increase in energy is likely to accelerate 
the dynamics, producing a negative change in C(0, t). 

The above equation shows that the three-point dy- 
namic susceptibility Xrit), up to a factor p/ (kBT^), mea- 
sures the spatial extent of the correlation between a fluc- 
tuation of energy at time and at some point 0, and the 
change in the dynamics occurring at another point r be- 
tween times and t. Clearly, XT{t) qualifies as a relevant 
quantitative indicator of space-time correlations in glass- 
formers, whose temperature (or density) dependence is 
now experimentally accessible. 

How does this alternative indicator relate to the dy- 
namic susceptibility X4(i)? On physical ground, it seems 



reasonable that energy fluctuations provide the main 
source of fluctuations in the dynamics, especially within 
a picture of slow activated relaxation with sizeable en- 
ergy barriers (recall that we consider the NVT ensem- 
ble). As a result, spatial correlations between changes in 
the dynamics at different points can be mediated by the 
energy fluctuations, as thoroughly discussed in Ref. [1^. 
This view is supported by an exact relation that can 
be derived between XT{t) and X4(i) [H, [H, [1^. From 
the Cauchy-Schwarz inequality, {SE{0)SC{0,t)) 



NVT 



{6E{0f)j^^^{SCi0,t)^)j^y^, and the thermodyn amic 
relation, N (^SE{0)'^) ^^^^ = ksT'^Cy, with c„ the con- 
stant volume speciflc heat per molecule in fc^ units, one 
indeed obtains 



ksT^ 

Cy 



[xr'^it)]' 



(8) 



so that the experimentally accessible response XT(t) can 
be used to give a lower bound on the well-studied suscep- 
tibility X4(^)- ^ rapidly growing energy-dynamics corre- 
lation as temperature decreases therefore directly implies 
a rapidly growing dynamics-dynamics correlation [33| . 

Actually, it can be shown that X4 and xr are more 
precisely related by 



xr^it) 



ksT'- 



\NVT 

Xt 



{t)\ 



xr^(t). 



(9) 



The dynamic susceptibility in the NVE ensemble, 
^NVEf^^y is a variance and therefore always positive, 
which of course is compatible with the inequality in 
Eq. ([HI). It is no more accessible to experimental mea- 
surements than Xi^'^i't)^ but it can nevertheless be com- 
puted in numerical simulations and in some theoretical 
approaches, as discussed in Refs. (33l. [34j. In particular, 
it was shown from simulations on two different models of 
glassforming liquids that for times of the order of (and 
in the range covered by the simulations) Xi^^ becomes 
small compared to (A:bT^/c«)(Xt^"^)^ when temperature 
is low enough and the dynamics becomes truly sluggish, 
which then makes the lower bound in Eq. ([H]) a good 
estimate of Xi^'^ ■ More importantly, the temperature 
dependences of Xi^^ and |xt^"^I ^o^' i ~ are found 
to be very similar, as can be proved exactly within the 
framework of the Mode-Coupling Theory. Furthermore, 
for Brownian dynamics Xi behaves as Xi^^ ■ These re- 
sults strongly support the claim that the dynamic spatial 
correlations in glassforming liquids are in fact fundamen- 
tally contained in the three-point dynamic susceptibility 
XT (see [13, [H, [13 and the discussion in Sec. Ill El be- 
low). In this respect, it is satisfying that XT{t = ^a) is 
found to be independent, for large Tq, of the type of dy- 
namics (Newtonian, Brownian or Monte-Carlo), as is the 
very phenomenon of glassy slowing down (38l . [Sa , [40l . [4]| . 
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C. Experimental conditions: NPT thermodynamic 
ensemble 

Most experimental situations in glassforming liquids 
and polymers correspond to the NPT ensemble. In order 
to carry out the analysis of the experimental data, one 
must first discuss the modifications brought about in the 
above presentation by switching the statistical ensemble 
from NVT to NPT. 

The counterparts of Eqs. ([6][9l) are simply obtained by 
replacing energy by enthalpy and constant volume spe- 
cific heat by constant pressure specific heat. For instance, 
one now has in place of Eq. ([7]) 

^T^'^'it) = d'rmO; 0)5c{r; 0, t))^pT, (10) 

where Xt^^(0 the derivative of F{t) with respect to T 
when TV and P are kept constant and where h{r, t) is the 
enthalpy density per molecule. Equation ^ is similarly 
changed into 

xr^ = ^(xrT+x4™, (11) 

Cp 

where the first term in the right-hand side is a lower 
bound for the nonlinear susceptibility Xa^'^ Xt^^ 
is, up to a factor p/ikBT"^), the integral of a correlation 
function between local fluctuations of the dynamics and 
of the enthalpy. 

Under constant pressure conditions, a change in tem- 
perature is accompanied by a change in density. As a re- 
sult, one can separate in the response of F{t) to a change 
of temperature at constant pressure the effect of temper- 
ature at constant density, described by Xt^"^ defined in 
Eq. ([HI), and the effect of density at constant temperature, 
characterized by x^^'^jt) = dF{t)/dp\T- This allows us 
to generalize Eq. (fTT|) to 

Xa — VXt ) +P'^B-LHT(Xp ) +Xa 

Cy 

where kt is the isothermal compressibility of the liquid. 
The last term of the right-hand side being positive, one 
can now use the sum of the two first terms as a lower 
bound to Xa^^ J which is slightly different from that con- 
sidered above, but which explicitly takes into account all 
sources of fluctuations. 

In any case, since both Xt ^^'^ X^^'^ can be experi- 
mentally obtained, the above expressions will allow us to 
evaluate the relative magnitudes of the contribution due 
to energy fluctuations and to density fluctuations in the 
space-time correlations, which represents an interesting 
issue per se. 

D. From dynamic susceptibilities to correlation 
volumes 

So far, we have presented a rigorous framework in 
which the spatially heterogeneous nature of the dynamics 



in glassforming materials is described through suscepti- 
bilities that involve spatial correlations associated with 
the dynamics. One expects that those correlations die 
out at large separation, so that their integral over space 
provides a measure of the characteristic "correlation vol- 
ume" at a given time t. It is then tempting, in light 
of Eqs. (0 El [TO]), to interpret an increase with decreas- 
ing temperature of cither x^i^) or kpT'^lxTit)] (for t of 
the order of Tq ) as a signature of a dynamic correlation 
lengthscale or correlation volume that grow as the glass 
transition is approached. 

Although the above considerations appear most plausi- 
ble and indeed, as will be discussed below, are supported 
by a series of numerical and theoretical arguments, it 
is also worth stressing what could go wrong in the cor- 
responding fine of thought. With X4(t) and xrit) one 
only has access to the integral over space of a correla- 
tion function, and two potential difficulties must be ad- 
dressed. First, the correlation function should go to zero 
at large distances. Although this may sound as a triv- 
ial requirement, it is more subtle than anticipated and 
it depends on the statistical ensemble under considera- 
tion [33| . It is a well-known feature, even in the simplest 
case of the static pair correlation [33| , that large distance 
limits of correlation functions in different ensembles dif- 
fer by terms of the order 1/A^ or 1/V. The difference is 
negligible in the thermodynamic limit, except when the 
function is inte gra ted over the whole volume, as it is in 
susceptibilities |42| . It turns out that this is not an issue 
when all the conserved quantities are let free to fluctuate. 
For one-component molecular liquids, this corresponds to 
the NPT ensemble in which both energy and density may 
fluctuate; on the other hand, the (negative) asymptotic 
background contribution has to be taken into account in 
the NVT ensemble, which we will therefore no further 
consider. 

Secondly, in the absence of more direct information 
on the spatial dependence of the correlation functions, 
extracting from their integral over space a "volume of 
correlation" (the derivation of a "correlation length" will 
be discussed below) requires an estimate of the ampli- 
tude of the correlation function to properly normalize the 
function. Indeed, an increase of X4 or kpT'^lxTl as T de- 
creases could be due not only to the extension in space of 
the associated correlations, the phenomenon one is look- 
ing for, but also to a variation of the typical amplitude 
of the correlations. Disentangling the two effects is thus 
a prerequisite. 

In the following, we focus on a timcscalc of the order of 
the a-rclaxation time , for which the susceptibilities X4 
and I XT I are found to reach a maximum. Relaxation on 
this timescale, at least at low temperature, is known to be 
due to rare dynamic events associated with a signiflcant 
apparent activation energy barrier (compared to ksT). 
This implies that the relevant dynamical events typically 
affect c(r, 0, t ~ r^) on the order of a fraction of the mean 
value of the correlation function. Therefore, provided the 
function C(0, t) is normalized by its mean value at < = 0, 
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F{Q) = (C(0, 0)), which we shall assume from now on, 
such changes should be of order one and should not vary 
much with temperature. 

The case of the dynamic susceptibility Xii^) is well doc- 
umented. One expects that the typical magnitude of the 
auto-correlation function {Sc{0; 0,t ^ Ta)6c{r; 0,t ^ Tq)) 
when r is in the vicinity of point is of order one when the 
dynamics is dominated by activated, intermittent pro- 
cesses. Hence Xii^ ^ Ta) is, up to a number of order 
unity, a measure of a correlation volume, or more prop- 
erly (see the density factor in Eq. ([5])) of a number of 
molecules that are dynamically correlated over a time 
span of the order of r^. More precisely, one can define 



^cori%4 = max{x4(i)}, 



(13) 



the maximum occurring for t ^ Tq,. 

Applying the same line of argument to the three-point 
susceptibility xt(^ ~ Tq) is a little more tricky because 
one now faces a cross-correlation function between en- 
thalpy fluctuation at time at a given point in space and 
fluctuation of the dynamics between t = and t ^ Tq in 
another place. 

First, the local enthalpy fluctuations that influence a 
dynamic event must have time to travel to the vicin- 
ity of the "active spot" . For instance, in the hydro- 
dynamic limit, long- wavelength and small- amplitude en- 
thalpy fluctuations follow a diffusive motion in a liquid, 
with a diffusion coefficient Dt\ the range of inffuence of 
such fluctuations must therefore be less than a length- 
scale of the order of V DtTo ■ However, because dra- 
matically increases with decreasing T whereas Dt has 
a much weaker dependence, this hydrodynamic length- 
scale is very large in viscous liquids and the upper bound 
it represents for dynamic correlations is of no real signif- 
icance. 

Secondly, the magnitude of the enthalpy fluctuations 
is temperature dependent, as shown by the thermody- 
namic relation N (^SH'^'^ = ksT'^Cp. Removing this un- 
wanted T-dependence requires to somehow normalize 
the variance of the enthalpy, e.g., to dcflnc a local di- 
mensionless fluctuation of enthalpy density, 5h[r), via 
N5H = py^ksT'^Cp J (PrSh{r). However, one may imag- 
ine that activated dynamic events associated with the a 
relaxation and leading to a change in Sc of order one 
involve enthalpy fluctuations more sizeable than aver- 
age, but rare. An appealing alternative is then to intro- 
duce another dimensionless enthalpy density fluctuation, 
6h{r), through 



(14) 



where Acp is the "configurational" part of the heat ca- 
pacity of a glassforming liquid, i.e., that in excess of 
the associated glass: this describes the contribution of 
those enthalpy fluctuations (indicated by the subscript a 
and whose typical size is given by N{6H\'^) = ksT^^Cp) 
that precisely disappear when the a relaxation is frozen 



out on the experimental timescale. Presumably, the re- 
maining fast components of the enthalpy fluctuations are 
then unable to induce the a relaxation, or at least are 
very weakly correlated with the a relaxation processes. 
(Nonetheless, and as usual, the temperature dependence 
of Acp can only be approximately determined experimen- 
tally by subtracting from the heat capacity of the liquid 
that of the crystal at the same temperature.) 

Note that this distinction between typical and rare 
fluctuations is especially important in the presence of lo- 
calized defects. Such defects, whose motion do not follow 
the standard hydrodynamic diffusion but is rather akin 
to that of "solitons" conserving their energy (or, at con- 
stant pressure, their enthalpy), have a negligible influence 
on the thermodynamics, hence on Cp(T), provided they 
arc dilute enough. In the case where such dilute defects 
are responsible for the slow relaxation of the system, it is 
clear that an estimate of the dynamic correlation volume 
based on characterizing the magnitude of the relevant en- 
thalpy fluctuations by means of Cp would be meaningless, 
whereas Acp in this instance would precisely capture the 
defect contribution and lead to a reasonable estimate, 
see Ref. [sj] for a discussion of defect models in this con- 
text. For glassforming liquids in general, the difference 
between Cp and Acp is not as dramatic, and in most cases 
remains of order of Cp itself [1^ . 

From the above discussion, we conclude that an es- 
timate of the number of molecules whose dynamics on 
the timescale of the a relaxation is correlated to a local 
enthalpy fluctuation is given by 




d^r {5h{O;0)5c{r;0,t ^ To,)) 



max{|xT(t)|}: 



(15) 



where Sh is the dimensionless fluctuation of energy den- 
sity introduced above and where we have used the fact 
that |xT(i)| is maximum for t ^ Ta [ssl . [33 |. This is 
the central formula which bolsters our following analysis 
of experimental data in terms of number of dynamically 
correlated molecules. 



E. Discussion 

Before moving on to the actual analysis of experimen- 
tal data, we would like to stress a number of points. 

(i) NcoriA versus A^corr.T A priori, iVcorr,4 and A^corr,T 

need not coincide, nor should they share the same tem- 
perature dependence. The former is the typical num- 
ber of molecules whose dynamics are correlated "among 
themselves", the latter the typical number of molecules 
whose dynamics are correlated "to a local fluctuation of 
enthalpy" . For Newtonian dynamics, and if the bound in 
Eq. ([8]) is nearly saturated, then these two numbers are 
actually quite different since: 



N,. 



(16) 
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This is indeed the result found within Mode-Couphng 
Theory near its dynamic singularity [s^l • This difference 
is however entirely due to the presence of one (or sev- 
eral) conserved variables. The primary mechanism cre- 
ating dynamic correlations is captured by Ncon-.T and is 
due to the fact that any local perturbation coupled to the 
slow dynamics (energy, density, composition, etc.) affects 
the dynamics over a certain volume, which grows as the 
glass transition is approached. This is measured by a cer- 
tain dynamic response Xc((|r— r'|,i) = {dc{r,t)/de{r' ,0)) 
[3^134. [ssj. Now, summing over r for a given point per- 
turbation at r' to define a correlation volume is identi- 
cal to fixing r but perturbing the system uniformly in 
space (while staying in a linear regime) [IHl . This is the 
fundamental reason why xt, simply obtained by taking 
the derivative of C with respect to temperature, is re- 
lated to a dynamical response. By the same token, the 
derivative of C with respect to any quantity coupled to 
dynamics is also a dynamic response; in the regime where 
these quantities become large, we expect that all of them 
(Xt, Xp: ■•■) should behave similarly as a function of time 
and temperature, but with different prefactors (see (ssj 
for a detailed proof of this point). We shall check this 
prediction on experimental data in Sec. V below. 

Following the above reasoning, the difference between 
-^corr.T and A^corr,4 comcs from the fact that if energy is 
conserved, a local energy fluctuation at r affects the dy- 
namics in the surrounding, but also propagates in space 
to create another correlated dynamical fluctuation else- 
where. This, in effect, transforms iVcorr.T hito N'^^rvT- 
This justifies that wc will mostly focus below on Ncor-c.T 
or Ncoi-T.p and less on Ncom-A- One should nevertheless 
keep in mind that some physical observables, such as the 
non-linear (cubic) response to an external field, are di- 
rectly related to iVcorr,4 [1^, and that both three-point 
and four-point susceptibilities are deeply connected [33 |. 

(ii) Prefactors — The uncertainty concerning the pref- 
actors involved in extracting correlation volumes or num- 
bers of correlated molecules, in particular the normaliza- 
tion problem discussed above, or the difference between 



F{t) then comes through that of the relaxation time Tq, 
F{t,T) = <^>{t/Ta{T)) with $(0) = 1, which leads in the 
NPT ensemble to 



Nr.. 




$'(1) 



d log Ta 



dlogT 



(17) 



where the prime on $ indicates a derivative. By using 
now Eq. p6|) and dropping the factor |$'(1)| that quanti- 
fies the stretching of the relaxation function and is thus of 
order one, one arrives at A'^corr^ {d log Ta / d logT)'^ , 

which is similar to Eq. ([5]). 

The phcnomenological derivation can thus be recast in 
a rigorous framework, which allows one to discuss the 
approximations made and, quite importantly, to under- 
stand the nature of the correlation that is probed. Essen- 
tially, the latter is the spatial correlation between local 
fiuctuations of the dynamics and of the enthalpy. Any 
interpretation in terms of Adam-Gibbs-like coopcrativ- 
ity regions [3| requires additional, speculative steps that 
seem of dubious validity, as discussed further in Sec. IVII 

(iv) Link with numerical simulations — A major advan- 
tage of having a rigorous formulation of the space-time 
correlations in liquids, relating dynamic susceptibilities 
to multi-point correlation functions, is that the heuris- 
tic considerations that we have developed in this and the 
previous subsections can be checked on models of glass- 
forming liquids. Computer simulations arc extremely 
useful in this context: one is guaranteed that the fea- 
tures of the dynamic heterogeneities that arc probed in 
experiments arc indeed the same as in simulations, ex- 
cept that the latter allow a much more detailed analysis, 
in particular concerning the full spatial dependence of 
the multi-point correlation functions, and not only their 
integral over space such as xt or Xi- The numerical in- 
vestigations so far carried out unambiguously confirm the 
interpretation of Eqs. [TB)) for estimating the spatial 
extent of the correlations in the dynamics of glassforming 
liquids 

(v) Lengthscales — The growth of the number of cor- 
using Cp and Acp, precludes a precise determination of related molecules, TVcon-.T, as one lowers the tcmpcra- 



the absolute value of those quantities. For instance, in 
the case of strong glassformers, the ratio Acp/cp may 
vary significantly, say from 1/2 to 1/10, between a fragile 
and a strong liquid. Choosing Cp or Acp, or even Acp/Cp 
as proposed by Donth [3^, then affects the relative esti- 
mates of the strong versus fragile glassformers. As found 
in Ref. [Ill, the order of magnitude of the resulting cor- 
relation volumes appears nonetheless physically sensible. 
Here, we shall focus on the variation with temperature 
which, as shown in the next section, is less sensitive to 
the estimate of the prefactors. 

(Hi) Equation 0^ recovered — The simple phcnomeno- 
logical analysis described in Sec. Ill Al ls recovered by as- 
suming that the shape of the relaxation function Fit) 
in the a-relaxation regime does not vary with tempera- 
ture, what is commonly referred to as "time-temperature 
superposition" . The only temperature dependence of 



ture, or equivalently of the associated correlation vol- 
ume obtained by dividing A^corr.T by the density, is a 
strong indication of a growing correlation length. How- 
ever, the precise relation is not straightforward. There 
is no reason indeed for the correlation volume(s) to be 
compact and therefore to scale as iVcon- ~ ■ Actually, 
if such a behavior were true for Ncoh.t, Eq. (|16p would 
lead to NcoYiA ~ C^! Extracting a Icngthscalc is thus a 
highly nontrivial task which requires additional informa- 
tion on the spatial dependence of the multi-point corre- 
lation functions involved in xt(0 and Xiit)- Again, such 
information is provided by numerical simulations and by 
model theoretical calculations (e.g., in the framework of 
the mode-coupling approach, where iVcorr.T ^ S,'^ in the a 
regime, or in kinetically constrained models). The main 
result of both theoretical models and numerical simula- 
tions is that a unique lengthscale appears to characterize 



8 



both XT{t) and X4:{t)j the associated correlation volumes 
being however noncompact [H, [13, ■ 



III. DATA TREATMENT AND CALIBRATION 

We have shown in the previous section that an esti- 
mate of the number of dynamically correlated molecules, 
NcoTr,T (and via the lower bound, A^corr,4), can be ob- 
tained by measuring the sensitivity of two-time correla- 
tors to temperature changes, see Eq. A standard 
tool to access such a two-time correlator in supercooled 
liquids is dielectric spectroscopy, and we have gathered 
published (and original) dielectric data on a range of 
glassforming liquids. In addition, in a few cases, we have 
also considered relaxation data obtained by other tech- 
niques, photon correlation spectroscopy, dynamic light 
scattering, optical Kerr effect, and neutron scattering. 
In this section, we review our operational protocol for 
treating the experimental data and we examine the un- 
certainties due to systematic or statistical errors. Read- 
ers interested only in the main physical results may skip 
this section and focus directly on the following ones. 




FIG. 1: (Color online) Three different estimates (i-iii) de- 
tailed in Sec. lnTBl for the fuU curve Tdx{<J,T)/dT for glyc- 
erol at temperatures T = 196.6, 201.6, 207.1, 211.1, 217.2, 
and 225.6 K (from left to right). Full lines are for method (i), 
dashed lines for method (ii), and points for method (iii). The 
three methods give consistent results, with small deviations 
that are well understood. 



A. Fitting the relaxation data 



B. Robustness of the temperature derivatives 



The quantity measured in dielectric spectroscopy is 
the linear susceptibility T), which we have described 
at each available temperature by an efficient em piri cal 
paramctrization, the Havriliak-Negami (HN) form |44l | or 
the one proposed by Blochowicz et al. [4^ on the basis of 
the extended generalized Gamma function (GGE). The 
temperature dependence of the HN or GGE adjustable 
parameters have subsequently been fitted to a high-order 
polynomial or to a Vogel-Fulcher-Tammann form (the ac- 
tual form of the function is not important as long as it fits 
well the data). From the HN or the GGE parametriza- 
tions, one can directly compute the associated relaxation 
function in the time domain [4^ and then obtain the 
three-point dynamic susceptibility XT{t) by simply tak- 
ing the derivative of the fitted temperature paramctriza- 
tion. 

In the case of the other experimental probes that we 
have considered, the data is directly provided in the time 
domain. For light and neutron scattering, we have fit- 
ted the relaxation function to the standard Kohlrausch 
stretched-exponential (KWW) function, whereas for the 
optical Kerr effect, we have used the paramctrization 
given in Ref. [i^. As for the dielectric data, we have 
interpolated the temperature dependence of the parame- 
ters by polynomial or Vogel-Fulchcr-Tammann functions 
and used the resulting formulas to compute the derivative 
with respect to temperature. 

In our operational protocol, we strictly limit ourselves 
to the domains of measurement, avoiding extrapolations 
as well as combinations of several sets of data. 



Provided one does not extrapolate out of the ex- 
perimentally accessed temperature domain, the above 
paramctrization scheme provides a very good descrip- 
tion of the temperature and time dependences of the 
measured relaxation functions. However, one must re- 
call that we are ultimately interested in a derivative with 
respect to temperature of a quantity whose characteristic 
timcscale varies by many orders of magnitude for small 
temperature changes. One must therefore study the ro- 
bustness of the outcome. 

As a first check, we have focused on the dielectric relax- 
ation of glycerol. We have measured the linear dielectric 
susceptibility of glycerol every 1 Kelvin for a set of tem- 
peratures above Tg « 190 K and below 232 K d^]. From 
the measured complex capacitance, C{tjj, T), we have ob- 
tained the following normalized dynamical quantity: 



X(^, T) = —7T-: = Re 

e(0) - Eoo 



C{uj) - C(oo) 
C(0) - C(cx)) 



(18) 



where where e'(a;) is the real part of the dielectric sus- 
ceptibility. Contrary to e(0), eoo has not been directly 
measured and has instead been obtained from a fit of 
C(a;,T) using the HN paramctrization. It turns out that 
eoo/e(0) is a small number whose temperature depen- 
dence is sufficiently weak to be irrelevant in the calcu- 
lation of the temperature derivative dx{oJ,T)/dT. We 
have computed the latter by three different methods. 
The two first ones are as described in Sec. IIII Al x(^) 
has been fitted to (i) the HN and (ii) the GGE forms, 
and the temperature dependence of the parameters has 



been fitted as well; dx{i^,T)/dT has then been com- 
puted from the parametrization of the T and uj depen- 
dences. In the procedure (ii), wc have directly taken the 
parameter values given by Blochowicz et al. [45| from ad- 
justment to their own dielectric measurements on glyc- 
erol. The third method (iii) consists in a direct calcu- 
lation of dx{(^,T)/dT through finite differences on the 
experimentally measured capacitance: dxii-o,T)/dT ~ 
[x(w, T+^)- x(w, T ~ ^)] /AT, with AT a small, 
finite temperature step. Note that we have considered 
the frequency domain instead of the time domain because 
this allows a direct use of the raw data when computing 
the finite differences, without further manipulations. 

The three different estimates obtained through (i-iii) 
are shown for several temperatures in Fig. [TJ The results 
are very similar. At a given T, 9x(w,T)/9T reaches a 
maximum, dx/dT\*, at a frequency to* which is close to 
the frequency uja ^ at which the imaginary part 
e"(a;) of the dielectric susceptibility is maximum. The 
values oj* and dx/dT\* are very close for the three meth- 
ods, with small differences which we now discuss. 

The comparison of methods (i) and (iii) reveals that 
performing finite differences of raw data does not change 
the value of lu*, but slightly underestimates the maxi- 
mum dx/dT\* by an amount which increases with AT. 
For AT = 1 K, the maximum value obtained through 
(iii) is lower than the two others by about 4%. The er- 
ror on dx(j^,T)/dT due to our experimental uncertain- 
ties on the dielectric measurement is estimated to vary 
from less than 1% at the maximum dx/dT\* to 5% when 
Tdx{i^,T)/dT is ten times smaller. In Fig. [1] it is worth 
noticing that method (iii) yields smooth curves despite 
the fact that no fitting procedure of the data is involved. 
This is consistent with errorbars of at most 1% at the 
frequency lu*. 

As for the comparison between methods (i) and (ii) 
the difference between the values of dx/dT\* is small, 
less than 3%, as shown in Fig. [T] The main difference 
lies in the value of uj*{T). Our values (following (i)) are 
smaller than those of Blochowicz et al. by typically 
20%. This amounts to decreasing all our temperature 
values by 0.6 K. Beyond a mere problem of thermometer 
calibration, this discrepancy may come from a difference 
in the purity of the glycerol samples, in particular their 
water content. It has been shown [4^ that a glycerol 
sample protected from any humidity yields, at a given 
T, values of Tq, larger than those of a sample "regularly 
exposed" to air humidity, corresponding approximately 
to a temperature shift of 1 K, at least in the limited 
temperature interval considered here. This is quite larger 
than the difference reported in Fig. [TJ To reinforce this 
interpretation, let us mention that in another series of 
experiments (not reported here) we have found values of 
uj*{T) closer to those of Ref. [45|. The difference between 
our two sets of data amounts to a 0.17 K shift, while the 
exact same thermometer was used in both experiments. 
The important point for our present purpose is that the 
difference in uj*(T) does not come with any significant 
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FIG. 2: (Color online) A^corr.T versus Tc/to on a log-log plot 
for glycerol (top) and o-terphenyl (bottom). Ncorr,T is calcu- 
lated by using Cp{Tg) (circles), Cp(T) (squares), Acp(rg) (dia- 
monds) and Acp{T) (triangles). The two upper (resp., lower) 
curves correspond to using Acp (resp., Cp); tq is arbitrarily 
set to 1 ps. 



difference in the maximum value of dx/dT\*. 

These detailed experiments performed on glycerol 
show that the temperature derivative involved in A'corr 
can be consistently obtained by using either finite dif- 
ferences (as also done in numerical work (33j). or a 
parametrization of the temperature and time (or fre- 
quency) dependences of the relaxation data through the 
HN or GGE descriptions. 



C. Dependence of Ncoii,t on thermodynamic inputs 



As discussed in Sec. |TT] and illustrated by Eq. (|15p . 
determining NcoTr,T (and consequently A'corr,4) requires 
thermodynamic input, primarily the isobaric specific 
heat Cp. For a few representative glassforming liquids for 
which enough data are available, we have investigated 
the two following points: (1) the quantitative difference 
in the estimate of A^coit.t and its temperature dependence 
which is introduced by using either Cp or Acp, the value 
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in excess to that of the associated crystal (see the dis- 
cussion in Sec. IIID|) : (2) the effect on the T-dependence 
of Ncoir.T that results from replacing Cp(T) (or Acp(T)) 
by a constant value, taken at Tg, Cp(Tg) (or Acp(Tg)). 
There are indeed a number of glassforming materials for 
which only the latter value is known. 

The results of our comparative study on two systems, 
glycerol [i^ and o-terphcnyl [s^, are displayed in Fig. [21 
where we plot Ncon-.T as a function of Tq/tq, where Tq. 
is the a-relaxation time extracted from the HN fit of 
the dielectric data and tq is an arbitrarily chosen micro- 
scopic time of 1 ps. One observes that the evolution is 
essentially the same in the different cases, using either 
Cp{T), Acp{T), Cp{Tg), or Acp{Tg). Replacing Acp{T) 
by a constant Acp{Tg) leads to a maximum deviation in 
Ncoir,TiT) of about 10—15% at high temperature around 
the melting point. Obviously, using Acp in place of Cp 
increases the absolute value of A'corr.T, by up to a fac- 
tor 2 in the case of o-terphenyl, but the whole shape of 
]^coit,t{T) is not significantly altered. Therefore, in what 
follows, we shall present results calculated with Acp(Tg), 
or, when not available, with Cp{Tg). 



D. Dependence of Ncoii,t on the experimental 
probe 

In our estimates of A'con- {T) we have focused on dielec- 
tric spectroscopy for which the largest data base is avail- 
able. One may however inquire how much would the 
results change if other experimental probes were used. 
On the one hand, one expects that spatial correlations of 
the dynamics arc an intrinsic property of the material, 
irrespective of the probe monitoring the dynamics. On 
the other hand, the estimates provided by Eq. ([TS]) and 
(fTB)) are not as straightforward as one would like (see the 
discussion in Sees. Ill Dl III E[) and they may depend on the 
probe chosen. This can be simply illustrated in the case 
where the relaxation functions obey "time-tcmpcrature 
superposition". As shown by Eq. (fT7)) . A^coir,T is pro- 
portional to I $'(1)1, which is essentially the stretching 
parameter P (it is equal to P/e if ^{x) is a stretched ex- 
ponential, ). If two probes are characterized by the 
same T-dependence of their relaxation time but by differ- 
ent stretc hing parameters, which is not uncommon (see 



e.g. Ref. |5l|), the resulting estimates of iVcorr.T then 



differ: the more stretched the relaxation (smaller /3), the 
smaller the value of A^coit.t- 

Matter can be worse if time-temperature superposi- 
tion does not hold. We illustrate this point by consider- 
ing several sets of relaxation data for liquid m-toluidine: 
low-temperature (near Tg) dielectric and photon cor- 
relation spectroscopy (PCS) [s^ measurements, high- 
temperature (above and around melting) neutron scat- 
tering data. Over the limited domains of temperature 
studied, both dielectric (s^ and neutron scattering relax- 
ation functions appear to satisfy time-temperature super- 
position with (accidentally) almost the same stretching 
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FIG. 3: (Color online) Ncoir,T versus Tg/T for m-toluidine 
evaluated from dielectric (triangles) and photon correlation 
(circles) spectroscopies close to Tg, and by coherent quasi- 
elastic neutron scattering (squares) at high temperature. (We 
have not considered the two lowest temperature data points 
for which the determination of the relaxation profile is too un- 
certain.) The inset shows the corresponding relaxation times 
(determined from the maximum of the imaginary part of the 
dielectric susceptibility, from a KWW fit of the correlation 
functions obtained from photon correlation spectroscopy and 
neutron scattering a.t q = 1.3 A^^). 



exponent (3 ~ 0.6 (the dielectric Kohlrausch stretching 
exponent is estimated from the fit in frequency domain 
according to [11]). On the other hand, the PCS data 
show a marked increase of the stretching, i. e. a decrease 
of /3, as Tg is approached [H^l- The resulting estimates of 
NcoTr,T are shown in Fig. [3] . 

The PCS data lead to a smaller value of Ncoii,t than 
the dielectric ones (by a factor of about 2) and to a slower 
growth as one lowers the temperature near Tg. This be- 
havior comes from the competition between the increase 
of log Ta and the decrease of /3 as T decreases. Simi- 
lar results could be expected for other molecular liquids 
due to analogous differences between dielectric and light- 
scattering dispersions [H, [sH] ■ 

The observation that average dynamic relaxations 
might depend on the probe is one of the puzzling ex- 
perimental features of the glass transition. One usually 
gets away by arguing that from different probes are 
nevertheless consistent with one another, while differ- 
ent shapes result from different microscopic observables 
(when they can be properly defined) . This becomes even 
more puzzling when these data are derived to obtain the 
value of iVcorr.T- If the observed averaged behaviors do 
not coincide, then it is no surprise that derivatives look 
even more different, as seen in Fig. [3l But this raises 
interesting questions if one wants to interpret the result 
in terms of a correlation volume: are there experimental 
probes that are "better" than other because more directly 



11 



coupled to molecular degrees of freedom? Can the high 
precision obtained in the dielectric spectroscopy experi- 
ments of Sec. IIIIBI be obtained using other probes? Is 
the discrepancy observed in Fig. [3] due to the compari- 
son of data with different resolutions? Can different mi- 
croscopic observables be correlated over different length- 
scales? Is it possible to extract a probe-independent cor- 
relation volume? More work is certainly needed to answer 
these questions. 



IV. TEMPERATURE DEPENDENCE OF Ncow 
FOR A RANGE OF GLASSFORMERS 

Having detailed how to estimate from experimen- 
tal data the temperature dependence of the number of 
molecules whose dynamics are spatially correlated and 
discussed the robustness of our method, we now present 
the results for a variety of glassforming liquids. The data 
shown in Figs. |4] and [5] represent the central outcome of 
our work. 

In Fig. m (top) we plot on a logarithmic scale Ncon-.T 
determined using Eq. (fTS)) versus Tc/tq, where is the 
a-relaxation time and tq is a microscopic time arbitrarily 
fixed to 1 ps (except for a couple of systems discussed be- 
low). Due to the absence of data for Acp for the last three 
systems, results for all substances in Fig. |4] are computed 
using Cp(Tg). The data include glycerol, o-terphenyl [s^ . 
m-toluidine |53 . salol [5^, propylene carbonate [4^, m- 
fluoroaniline |45|. propylene glycol 45 1, and decaline [57| 
(dielectric measurements), as well as B2O3 153 fphoton 
correlation spectroscopy), BPM [l^ and salol |59l| (opti- 
cal Kerr effect), and m-toluidine (neutron scattering). 
We have added for comparison simulation results ob- 
tained on two very different model glassformers, namely 
the BKS model for silica [11] and a binary Lennard-Joncs 
mixture [1^. In those two cases, the correlation functions 
used to probe the dynamics are self-intermediate scat- 
tering functions, and temperature derivatives are mea- 
sured using finite differences of measurements performed 
at nearby temperatures, T and T -I- ST, with ST small 
enough that linear response applies. For BKS silica, we 
take To = 1 ps, while we take the standard Lennard-Joncs 
time unit as a microscopic timescale, tq = ^ma^/{A8e), 
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where m, a and e respectively represent the mass of the 
particles, their diameter, and the depth of the Lennard- 
Jones potential; using argon units one gets tq « 0.3 ps. 
Finally, we have included a colloidal hard-sphere system, 
simply reporting the results already published in Ref. (3^ 
and choosing tq = 1 ms. 

A key feature in the results presented in Fig. 2] is 
that Ncorr,T (and via Eq. (jTB]), A^corr,4(T)) does increase 
as temperature decreases or the relaxation timescale in- 
creases. The extent of spatial correlation in the dynam- 
ics therefore grows as one approaches the glass transi- 
tion [3^ . 

Although there is an unfortunate shortage of data cov- 
ering the whole temperature range from above melting 
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FIG. 4: (Color online) Top; Evolution of Ncott,t for the 
different materials indicated when the glass transition is ap- 
proached. The results are shown on a logarithmic scale as a 
function of Tc/tq, where tq is a microscopic timescale (see the 
text). The full line is from Eq. p9|l and describes a crossover 
from a power law scaling at high temperatures (small r^) to 
a logarithmic growth close to the glass transition, Eq. p9|l . 
Bottom: Using the freedom left by unknown normalizations 
of order unity we obtain a better collapse of the data onto the 
fit. 



down to Tg, Fig. |4] clearly indicates that Ncoiy.t grows 
faster when Tq is not very large. The initial, fast growth 
observed in A^corr occurs at high temperature, close to the 
onset of slow dynamics. A power law fit, A^corr,T ~ Tq^''', 
in this regime {tc/tq < 10^) leads to values 7 = 2 — 3. 
Interestingly, these are consistent with the predictions 
of Mode Coupling Theory, see [H, [3^ . The growth be- 
comes much slower as Tg is approached. More quantita- 
tively, if one wants to fit with a (second) power law in 

the regime Tq/tq e [10^, 10^^], one finds Ncoii,t ^ , 
with 7' > 20 and small variations among the different liq- 
uids. This means that when the dynamics slows down by 
6 decades in time, A'coit.t merely increases by a factor 
about 2. This slow growth is actually broadly consis- 
tent with activation-based theoretical approaches which 
all predict some sort of logarithmically slow growth of 
dynamic lengthscales [60I . [6ll . [g^ . [63| . In order to cap- 
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FIG. 5: (Color online) NcorT,T/Ncoy:y:,T{Tg) (top) and Tc/tq 
(bottom) versus Tg/T for different glass-forming liquids (in- 
dicated in the Figure) on a logarithmic scale; Tg is defined by 
Ta,{Tg) = 10 s. 



ture the idea of a fast initial growth of A^corr,T when 
is small, followed by a logarithmically slow growth, we 
have empirically fitted the data to the formula: 



A 



N. 



corr,T 







exp 



iV, 



corr,T 

No 



(19) 



and found rough agreement (see the full line in Fig. 3]) 
with A = A, No = 0.8, 7 = 2, and V = 1-4. (Note 
that very similar ideas have been used in spin glass stud- 
ies (63|.) These values are unfortunately only indicative 
because a set of rather different values allows one to fit 
almost equally well the curves. It would be interesting to 
be able to narrow down the range of acceptable values for 
ip in order to test more stringently the various activation- 
based theories found in the literature [13, HU, IH, . 

A second interesting feature of the results presented in 
Fig. m is that the data for different materials fall remark- 
ably close to one another, despite the various caveats 
mentioned above about prefactors, normalization etc.. 



glassformer 


NcOYY,T{Tg) 


A'"corr,4(T'9) 


Acp 


Cp 


glycerol 


8.2 


67.2 


90.5 


175 


o-terphenyl 


9.2 


84.6 


112.7 


333.7 


salol 


8.6 


74.0 


118 


320 


propylene carbonate 


15.9 


252.8 


75.4 


164 


m-fluoroaniline 


12.9 


166.4 


86 


161 


propylene glycol 


9.4 


88.4 


67.2 


150 


B2O3 


10.1 


102.0 


40 


131 


m-toluidine 


11.7 


136.9 


89.5 


192.6 


decaline 


11.8 


139.2 


64 


133 



TABLE I: Values of Ncoii,T{Tg) for various glassforming liq- 
uids, and corresponding lower bound on Acorr,4(Tg). The val- 
ues are computed with Acp{Tg). For completeness we report 

-1. mol-i. 



which might affect the vertical axis in this figure. A 
signature of such caveats is the fact that we find val- 
ues of iVcorr,T at high temperature which fall below 1. 
Investigating in detail the reasons of this behavior and 
improving the normalization procedure, however, would 
require the knowledge of the spatial distribution of the 
dynamical correlations, which is for the moment out of 
reach experimentally. Moreover our choice of a micro- 
scopic timescale (mainly 1 ps) is also somewhat arbi- 
trary. Clearly, a better collapse of the different curves 
could be obtained by making use of the freedom offered 
by these unknown normalizations. In Fig. 3] (bottom) 
we rescale vertical and horizontal axis by factors or or- 
der unity (i.e. we only allow data shifts within a decade, 
very often much less) in order to obtain a better collapse 
along the fit described above. We find that all data col- 
lapse rather well, suggesting that Eq. captures the 
physics of glassformers very well. 

We show in Fig. O (top) the same data as in Fig.[51 nor- 
malized by the value of iVcorr,T at Tg, versus Tg/T (ex- 
cluding the simulation results and the colloidal system, 
for which Tg cannot be determined) . Table [T] shows the 
different values of A^corr.TCTg). Recall that estimates for 
A^corr,4(T) are simply obtained through Eq. which 
on a logarithmic scale merely amounts to rescaling the 
y-axis by a factor of 2. 

In Fig. [5] we also compare the celebrated Angell plot 
(bottom) that illustrates the relative fragility of glass- 
formers HI] to the relative change of Ncott,t{T) (top). 
We do not find any systematic correlation between the 
two plots, i.e., between the rapidity of the increase of 
Ncorr,T{T) and the fragility (with the notable exception 
of decaline which is also the most fragile material con- 
sidered). Indeed, the results for all the systems show 
a rather similar behavior, characterized above all by a 
fairly modest growth of A^corr,T- We note in passing that 
since Acp is taken here as constant and replaced by its 
value at Tg, Eq. (fT7|) would on the contrary predict that 
NcoiiiT) follows the fragility pattern. The fact that this 
is not what is observed is mainly due to the breakdown 
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of the time-tempcraturc superposition property and the 
resulting effect of the shape (stretching) factor. 

Finally, and keeping in mind the theoretical and prac- 
tical uncertainties on the determination of the absolute 
value of A'corr, we note that the values of Ncorr, T{Tg) 
for all studied glassformers are similar, all between 8.2 
and 15.9 as shown in Table U Contrary to what has 
been done in Ref. [3^ . we have not attempted here to 
convert this number into a precise length scale (see the 
discussion above), nor to "normalize" Ncori-,TiTg) (or 
Ncon-.iiTg) ~ A^coir,T(T'g)^) by usiug an effective num- 
ber of "beads" per molecule to account for the difference 
in complexity of the various molecules. The relative inde- 
pendence of A^corr(Jg) upon fragility is broadly consistent 
with Wolynes's version of the Random First Order The- 
ory [b^, [6y] , although the relation between iVcon- and the 
typical size of the mosaic state is far from obvious - see 
the conclusion for a deeper discussion of this point. 

V. RELATIVE EFFECTS OF TEMPERATURE 
AND DENSITY FLUCTUATIONS ON 
DYNAMICAL CORRELATIONS 

The last question we want to address is the relative 
effects of temperature and density fluctuations on dy- 
namical correlations. In Sec. Ill C[ we have shown that 
the relative influence of temperature and density fluctu- 
ations on Xi is represented respectively by the contribu- 
tions ksT^/c^ixT^^? and p'^kBTKriXp^^? ■ We have 
also argued in Sec. Ill El that both terms in fact represent 
a (squared) dynamic response and should behave simi- 
larly as a function of temperature, provided both T and 
p locally couple to the dynamics. It is therefore of in- 
terest to compare these two contributions, which we now 
do for three molecular liquids characterized by different 
fragilities: glycerol, m-toluidinc, o-tcrphcnyl. 

We introduce a number of simplifications in our pro- 
tocol. (1) We assume that the property of "time- 
temperature superposition" holds for the temperature 
range under study. This allows us to write the relax- 
ation function as F{t,T) = <^{t/Ta{T)) with $(0) = 1 
and no explicit T-dependence elsewhere than in the re- 
laxation time Ta{T). (2) We consider the thermodynamic 
coefficients Cp, Ct,, kt, not the excess with respect to the 
crystal values, which usually are not all accessible. The 
data for Cp{T) and the e quat ion of state are taken from 
the literature @, M MM M M 

and Cv is obtained 

from the relation c„ = c„ -, where a„ is the iso- 

baric coefficient of expansion. (3) We use the scaling 
form onto which the density and the temperature de- 
pendences of the a-rclaxation time and the viscosity of 
many glassforming liquids and polymers have been found 
to collapse [tiI. [7a|. 

T„(p,T) = /[e(p)/fcBT], (20) 

where f{x) is a species dependent, scaling function and 
e(p) is a density-dependent (activation) energy scale. 
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FIG. 6: (Color online) A'^corr,4 versus Tg/T for glycerol, m- 
toluidine, and o-terphenyl (top to bottom) calculated via the 
lower bounds given in eg. ill! (squares) and 112! f diamonds) . 
Also displayed is the contribution due to energy fluctuations 
(circles). The inset shows the ratio R versus Tg/T (see eq. l22p 
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With the above simplifications, and dropping the irre- 
ducible contribution coming from Xi^^ i one finds from 
Eq. ini): 

iVeo„..4(p, T) ~ ^'{If^^ (1 + pTnw^) , (21) 

Cv 

where ra^ = — ^logrQ/^logTlv- is the "isochoric 
fragility" , which depends on p and T only through Tq (71| 
and = 91oge(p)/91ogp is found roughly constant 
over the range of densities involved at atmospheric pres- 
sure [tH . [t^ . The ratio of the contribution of the density 
fluctuations over that of the energy fluctuations is then 
given by: 

R = pTKTC„x^. (22) 

As shown in the insets of Fig. [6l we find that R is roughly 
constant and close to 0.4 for glycerol, varies as T de- 
creases from 0.7 to 0.6 for m-toluidine (note that in this 
case only values close to Tg are reported), and from 3 
to 2 for o-terphenyl. These results are compatible with 
our expectation, detailed in Sec. Ill El that density and 
temperature contributions mirror the same physical phe- 
nomenon, namely the anomalous susceptibility of the dy- 
namics to a local perturbation, and hence should share 
a similar temperature dependence. Thus, although the 
temperature contribution is dominant the other contribu- 
tion behaves similarly as the temperature is lowered. In 
particular the ratio of the two terms changes modestly 
compared to their overall variation, except perhaps for 
o-tcrphenyl. It is known that in the latter case, density 
plays a more important role, as noticed for the slowing 
down of the relaxation [t^I- It is however much less so 
in the other liquids where the contribution of the energy 
fluctuations dominate, especially at low temperature. 

As a byproduct of the present analysis, we can also 
evaluate the relative quality of the two lower bounds of 
Xi^'^ that are provided either by Eq. ([TTjl with Xi^^ 
neglected or Eq. with Xa^^ neglected. The result 
for the three liquids already considered above is shown 
in Fig. [S] Because it takes more properly into account 
the relevant sources of dynamic fluctuations, the lower 
bound calculated with Eq. (fT2|) is better than that ob- 
tained from Eq. pT|) . Although we do not have ac- 
cess to Xi^^ and Xa^^ individually, we can estimate 
their difference (evaluated as before for the time t ^ Ta 
at which Xa^'^ is maximum). At Tg, we find that 
ixr" - xr'')/{xr^ - xr") is 0.15, 0.06, 0.31 for 
glycerol, m-toluidine, and o-tcrphenyl, respectively. 

VI. CONCLUDING REMARKS 

Our central result has been to provide direct evidence 
from experimental data for an increase of spatial corre- 
lations in the dynamics of glassforming liquids as one 
lowers the temperature [3^ . Although the precise rela- 
tion to a lengthscale is not straightforward (as discussed 



in Sec. Ill Pp . this nonetheless shows the existence of a 
length that grows as the glass transition is approached. 
Moreover, our results establish an important character- 
istic of this phenomenon: despite the dramatic slowing 
down of the relaxation in supercooled liquids, the growth 
of Neon, and of the associated lengthscale, is fairly mod- 
est, especially below the onset of slow dynamics. This 
is compatible with an activation-based picture of the dy- 
namics in supercooled liquids. At higher temperature, 
however, the growth of this lengthscale appears to be 
faster, and not inconsistent with the predictions of Mode- 
Coupling Theory. 

The modest size reached by the correlation volume at 
Tg also conveys a two-fold message for future studies on 
the glass transition. (1) Experiments will presumably 
never be able to probe lengthscales as large as the ones 
measured for standard critical phenomena (where time 
and length scales are generically related in a power law 
fashion [23 )• As a consequence, it is unlikely that the 
growing of the lengthscale may provide indisputable evi- 
dence that the slowing down is driven by an underlying 
phase transition, whatever its nature. Similarly, theoret- 
ical results that neglect subleading contributions to the 
asymptotic critical behavior will have to be taken with 
a grain of salt. (2) Despite the huge gap in timescales 
between experimental and numerical work, the probed 
lengthscales are not very different, which indicates the 
ability of numerical simulations to capture some ma- 
jor aspects of the slow dynamics in glassforming mate- 
rials. A similar point was made in the context of spin- 
glasses [H,!?!- 

We have already stressed in this article the direct con- 
nection between A^corr and the dynamic heterogeneities 
observed in glassforming systems: A'corr can be taken as 
the number of molecules participating in a typical dy- 
namic heterogeneity. One may however wonder if Neon 
and the associated lengthscale relate to other aspects of 
the slowing down of relaxation leading to the glass tran- 
sition. Among these aspects are usually mentioned the 
"cooperative" nature of the a-relaxation 0, d, 01 and the 
stretching of the relaxation functions, often associated 
with a spread of local relaxation times [1, 0j Hi • 

A connection between TVcorr and what is usually im- 
plied by "cooperativity" seems to us dubious, and cer- 
tainly more work is needed to clarify this issue. Co- 
operativity, in the context of thermal activation, means 
that the effective activation energy is the sum of the el- 
ementary barriers associated with the motion of the ob- 
jects that evolve cooperatively. In the Adam-Gibbs pic- 
ture [1] these objects form a compact cluster, whereas in 
the frustration-limited domain picture [6l| and presum- 
ably also in the random first order transition approach 
near the thermodynamic singularity Tr- [60|.[63| these ob- 
jects form (possibly fractal) domain walls. In all cases, 
however, the total effective activation energy is directly 
related to the number of correlated molecules, say A^coop, 
through 

A(T) = ^eoop(r)''Ao, (23) 
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where Aq represents some elementary barrier that may 
or may not be temperature dependent, and is a certain 
exponent accounting for possible renormalization effects. 
As a result, t„ oc exp{N^^^^Ao/{kBT)). Now, using the 
above definition of A^con-.T in Eq. (jl5p and considering 
for illustration the simple picture in which 9 = 1 and 
Ao is independent of temperature (as in the traditional 
Adam-Gibbs approach), one obtains 



Nr. 



Ao 



Acp ksT """^ 



aiogA^e, 



aiogT 



(24) 



Since Acp is weakly temperature dependent, this expres- 
sion implies that A'corr (A'corr.T and, via Eq. lfT6|) A'corr,4 
as well) increases as temperature decreases even when 
A^coop is temperature independent (and equal to one, for 
that matter, which is tantamount to an absence of "coop- 
erativity"). This discussion underlines that, contrary to 
the concept of "dynamic heterogeneity" that can be prop- 
erly quantified, the notion of "cooperativity" in glass- 
forming liquids lacks an operational definition that could 
then lead to concrete measurements of Ncoo-p- 

We also point out that no meaningful relation seems to 
emerge between Ncotr.T and the stretching of the relax- 
ation. Here, we do not refer to the way the stretching pa- 
rameter enters in some expressions of A"corr,T and to the 
fairly unintuitive behavior that follows (see Sec. IIII D[) . 
We rather consider the common view of the stretching as 
a phenomenon which is attributed to the spatial hetero- 
geneity of the dynamics (and initially triggered the whole 
body of work on dynamic heterogeneities) and which is 
connected to the "fragility" of the glassforming materi- 
als [t^. The main argument supporting the absence of 
link is the example of a simple system in which relax- 
ation decays exponentially in time with a characteris- 
tic time that follows an Arrhcnius temperature depen- 
dence. There is no stretching (and no "fragility"), and 
yet, A^corr.T increases as temperature decreases. It should 
be stressed again that our conclusions concerning Ncott,t 
apply to Acorr.4 as well, since the two are intimately con- 



nected. 

An additional illustration of the above points is pro- 
vided is the case of the Fredrickson- Andersen (FA) mod- 
els [t^. fa models can be seen as defect models, where 
stochastic rules for the dynamical evolution of the defects 
are postulated. In the simplest form of the model (the 
"one-spin facilitated FA model"), the relaxation of the 
system proceeds by the thermally activated diffusion of 
point defects throughout the system. The average relax- 
ation time in the system behaves in an Arrhenius manner, 
while two-time correlation functions decay exponentially, 
at least in three spatial dimensions. Such a local relax- 
ation process can hardly be said to be "cooperative" . Yet, 
dynamics is very heterogeneous, in the sense that local 
dynamics is correlated on distances of the order of the 
typical distance between point defects, which can indeed 
become very large at low temperature. In this case, X4(0i 
XT, and A^corr all grow very strongly as temperature is 
decreased 0, [zllzl ■ 

For all the above reasons, we conclude on a caution- 
ary note: one should be careful in comparing numerical 
and experimental values of A^corr to theoretical predic- 
tions made by theories which do not directly investigate 
multi-point susceptibilities. This remark does not ap- 
ply to mode-coupling theory and kinetically constrained 
models, but does apply to most activation-based theories, 
such as the the Adam-Gibbs description [3], random-first 
order transition theory [b^, [1^ , or the frustration-limited 
domain approach (6l| . 
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